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Abstract. Given a simple polygon P in the plane, we present new al- 
gorithms and data structures for computing the weak visibility polygon 
from any query line segment in P. We build a data structure in 0(n) 
time and 0(n) space that can compute the visibility polygon for any 
query line segment s in O(fclogn) time, where k is the size of the visibil- 
ity polygon of s and n is the number of vertices of P. Alternatively, we 
build a data structure in 0(n 3 ) time and 0(n 3 ) space that can compute 
the visibility polygon for any query line segment in 0(k + logn) time. 



1 Introduction 

Given a simple polygon V of n vertices in the plane, two points in V are visible 
to each other if the line segment joining them lies in V . For a line segment s in 
V, a point p is weakly visible (or visible for short) to s if s has at least one point 
that is visible to p. The weak visibility polygon (or visibility polygon for short) of 
s, denoted by Vis(s), is the set of all points in V that are visible to s. The weak 
visibility query problem is to build a data structure for V such that Vis(s) can 
be computed efficiently for any query line segment s in V . 

This problem has been studied before. Bose et al. [2] built a data structure 
of 0(n 3 ) size in 0(n 3 log n) time that can compute Vis(s) in 0(k log n) time for 
any query, where k is the size of Vis(s). Throughout this paper, we always let k 
denote the size of Vis(s) for any query line segment s. Bygi and Ghodsi [3l gave 
an improved data structure with the same size and preprocessing time as that 
in [5] but its query time is 0(k + logn). Aronov et al. [1] proposed a smaller 
data structure of 0(n 2 ) size with 0(n 2 logn) preprocessing time and 0(fclog 2 n) 
query time. Table [T] gives a summary. If the problem is to compute Vis(s) for a 
single segment (not queries), then there is an 0(n) time algorithm 
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Table 1. A summary of the data structures. The value k is the size of Vis(s) for any 
query segment s. 



Data Structure 


Preprocessing 


; Time Size 


Query Time 


m 


0(n 3 log n) 


0(n 3 ) 


0(fc log n) 


3 


0(n 3 log n) 


0(n 3 ) 


0(fc + logn) 


m 


0(n 2 log n) 


0(n 2 ) 


0(felog 2 n) 


Our Result 1 


0(n) 


O(n) 


O(felogn) 


Our Result 2 


0(n 3 ) 


0(n 3 ) 


0(k + logn) 



1.1 Our Contributions 

In this paper, we present two new data structures whose performances are also 
given in Table [1] Our first data structure, which is built in 0(n) time and 
O(n) space, can compute Vis(s) in O(fclogn) time for any query segment s. 
Comparing with the data structure in pQ, our data structure reduces the query 
time by a logarithmic factor and uses much less preprocessing time and space. 

The preprocessing time and size of our second data structure are both 0(n 3 ), 
and each query takes 0(k + logn) time. Comparing with the result in [3], our 
data structure has less preprocessing time. In addition, our solution, which is 
based on the approach in [2], is much simpler than that in [3]. Further, our 
techniques explore many geometric observations on the problem that may be 
useful elsewhere. For example, we prove a tight combinatorial bound for the 
"zone" in a line segment arrangement contained in a simple polygon, as follows, 
which is interesting in its own right. 

Let S be a set of line segments in a simple polygon V such that both end- 
points of each segment of S are on dV (i.e., the boundary of V). Let A be the 
arrangement formed by the segments in S and the edges of &P. For any line seg- 
ment s in V (the cndpoints of s need not be on dV), the zone of s, denoted by 
Z(s), is defined to be the set of all faces of A that s intersects. For each edge of 
any face in A, it either lies on a segment of S or lies on dV. Let A be the number 
of edges of the faces in Z(s) each of which lies on a segment of S. We want to 
find a good upper bound for A. By using the zone theorem for the general line 
segment arrangement jS], we can easily obtain A — 0(\S\a(\S\)), where a(-) is 
the functional inverse of Ackermann's function |13j . In this paper, we prove a 
tight bound A — 0(m), where m < \S\ is the number of segments in S each of 
which contains at least one edge of the faces in Z(s). An immediate application 
of this result is that we obtain an efficient query algorithm for our second data 
structure. Since combinatorial bounds on arrangements are fundamental, this 
result may find other applications as well. 

The rest of this paper is organized as follows. In Section^ we review some 
geometric structures and a query algorithmic scheme that will be used by the 
query algorithms of both our data structures. We will also give a "ray-rotating" 
data structure in Section^ which is needed by our first data structure in Section 
[3] In Sections |3] and Q] we present our first and second data structures, respec- 
tively. As a by-product of our second data structure, the combinatorial bound 
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Fig. 1. Illustrating a window vw of p. 



Fig. 2. Illustrating the two critical con- 
straints vp^ and upZ defined by the two 
mutually visible vertices u and v. 



of the zone mentioned above is also given in Section 2) Section [5] concludes the 
paper. 

2 Preliminaries 

In this section, we review some geometric structures and discuss an algorithmic 
scheme that will be used by the query algorithms of both our data structures 
given in Sections |3] and 0J We will also give a "ray-rotating" data structure in 
Section 12.11 which is needed by our first data structure in Section [3] 

For simplicity of discussion, we assume no three vertices of V are collinear; 
we also assume for any query segment s, s is not collinear with any vertex of V 
and each endpoint of s is not collinear with any two vertices of V. As in [l|2j . 
our approaches can be easily extended to the general situation. 

Denote by dV the boundary of V. The visibility graph of V is a graph whose 
vertex set consists of all vertices of V and whose edge set consists of edges 
defined by all visible pairs of vertices of V. Here, two adjacent vertices on dV 
are considered visible to each other. In this paper, we always use K to denote 
the size of the visibility graph of V . Note that K — 0(n 2 ) and K — Q(n). The 
visibility graph can be computed in O(K) time [T4j . 

We introduce the visibility decomposition of V |l|2j . Consider a point p in 

V and a vertex v of V. Suppose the line segment pv is in V, i.e., p is visible 
to v. We extend pv along the direction from p to v and suppose we stay inside 

V (when this happens, v must be a reflex vertex). Let w be the point on the 
boundary of V that is hit first by our above extension of pv (e.g., see Fig.Q]). We 
call the line segment vw the window of p. The point p is called the defining point 
of the window and the vertex v is called the anchor vertex of the window. It is 
well known that the boundary of the visibility polygon of the point p consists of 
parts of dV and the windows of p |1I2| . If the point p is a vertex of V, then the 
window vw is called a critical constraint of V and p is called the defining vertex 
of the critical constraint. For example, in Fig. [2j the two critical constraints wpZ, 
and vp v are both defined by u and v; for up u , its anchor vertex is u and its 
defining vertex is v, and for vp^, its anchor vertex is v and defining vertex is u. 
It is easy to see that the total number of critical constraints is 0(K) because 
each critical constraint corresponds to a visible vertex pair of V and a visible 
vertex pair corresponds to at most two critical constraints. 
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As in |1I2) . we represent the visibility polygon Vis(s) of a segment s by a 
cyclic list of the vertices and edges of V in the order in which they appear on 
the boundary of Vis(s), and we call such a list the combinatorial representation 
of Vis(s) pp. With the combinatorial representation, Vis(s) can be explicitly 
determined in linear time in terms of the size of Vis(s). Our query algorithms 
given later always report the combinatorial representation of Vis(s). 

The critical constraints of V partition V into cells, called the visibility decom- 
position oiV and denoted by VT>("P). The visibility decomposition VD(V) has a 
property that for any two points p and q in the same cell of VD(V), the two vis- 
ibility polygons Vis(p) and Vis(q) have the same combinatorial representation. 
Also, the combinatorial representations of the visibility polygons of two adjacent 
cells in VD(V) have only O(l) differences. The visibility decomposition has been 
used for computing visibility polygons of query points (not line segments) [1I2| . 

Consider a query segment s in V . In the following, we discuss an algorithmic 
scheme for computing Vis(s). Denote by a and b the two endpoints of s. Suppose 
we move a point p on s from a to b. We want to capture the combinatorial 
representation changes of Vis(p) of the point p during its movement. Initially, p is 
at a and we have Vis(p) = Vis(a). As p moves, the combinatorial representation 
of Vis(p) changes if and only if p crosses a critical constraint of V [112] . Vis(s) 
is the union of all such visibility polygons as p moves from a to b. Therefore, to 
compute Vis(s), as in |H2j . we use the following approach. Initially, let Vis(s) = 
Vis(p) — Vis(a). As p moves from a to b, when p crosses a critical constraint, 
either p sees one more vertex/edge, or p sees one less vertex/edge. If p sees one 
more vertex/edge, then we update Vis(s) in constant time by inserting the new 
vertex/edge to the appropriate position of the combinatorial representation of 
Vis(s). Otherwise, we do nothing (because even though a vertex/edge is not 
visible to p any more, it is visible to s and thus should be kept; refer to [2] for 
details). 

The above algorithm has two remaining issues. The first one is how to com- 
pute Vis(a) of the point a. The second issue is how to determine the next critical 
constraint that will be crossed by p as p moves. Each of our two data structures 
given in Sections [3] and |4] does some preprocessing such that the corresponding 
query algorithm can resolve the above two issues efficiently 

2.1 The Ray-Rotating Queries 

Our first data structure in Section [3] needs the following ray-rotating queries. 
Given any ray p whose origin z is in V, the ray-rotating query asks for the first 
vertex of V visible to z that will be hit by p when we rotate p clockwise (or 
counterclockwise) around z (e.g., see Fig. [3]). By making use of the ray-shooting 
data structures [51611 1115] and the two-point shortest path query data structure 
[12] , we obtain the following result. 

Lemma 1. A data structure can be built in 0(n) time and 0(n) space such that 
each ray-rotating query can be answered in O(logn) time. 
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Fig. 3. Illustrating the ray-rotating query Fig. 4. Illustrating the proof for Lemma[T] 

for p: The vertex v is the first visible vertex The two dotted paths are shortest paths 

to z that will be hit by p if we rotate p from z to ui and Ui, respectively, 
clockwise around z. 



Proof: Consider any ray p whose origin z is in V. Without loss of generality, 
assume p is horizontally rightwards. Let v* be the sought vertex for the ray- 
rotating query on p, i.e., v* is the first vertex of V visible to z that will be hit 
by p when we rotate p clockwise around z (the case of counterclockwise rotation 
can be done similarly). 

In the preprocessing, we compute a ray-shooting data structure [516111115] 
in 0(n) time and space, such that given any ray with the origin in V, the first 
point on the boundary of V hit by the ray can be found in O(logn) time. We 
also compute a two-point shortest path query data structure [12] in 0(n) time 
and space, such that given any two points p and q in V, the shortest path length 
between p and q can be computed in O(logn) time and the path itself can be 
found in additional time proportional to the number of turns along it. 

Our query algorithm for finding v* works as follows. 

First, we use the ray-shooting data structure to find in 0(log n) time the first 
point q* on the boundary of V hit by p; the edge of V containing q* is also known 
immediately from the ray-shooting query. If q* is a vertex of V, then v* = q* 
and we are done; otherwise, let the end vertices of the edge of V containing q* 
be iti and ui (e.g., see Fig. [4j. Let %% be the shortest path in V from z to u%, 
and similarly, let -k-i be the shortest path from z to ui- Since z is visible to q* 
on U1W2, the region bounded by tt\ , n2 , and u\U2 is a funnel [11) 12) 1 7j , with z 
as the apex (e.g., see Fig. |4j. Recall that p is horizontally rightwards; one vertex 
of iti and U2 must be below the line containing p. Without loss of generality, let 
Ui be below the line containing p. Let v be the vertex on 7Ti that is connected 
to z by a line segment on m, i.e., ~zv is the first edge of m (e.g., see Fig. [4]). 
Note that v = iti is possible, in which case 7Ti is the line segment zu\. An easy 
observation is that the sought vertex v* is exactly the vertex v. By using the 
two-point shortest path data structure |12j on z and u\, the vertex v can be 
easily found in O(logn) time because ~zv is the first edge of 7Ti. 

Therefore, the sought vertex v* for the ray-rotating query on p can be found 
in 0(log n) time. The lemma thus follows. □ 
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Fig. 5. Illustrating the principle child w of v in T(p). 

3 The First Data Structure 

Our goal is to compute Vis(s) for any query segment s. Again, let s = ab. As dis- 
cussed before, we need to resolve two issues. The first issue is to compute Vis(a). 
For this, as discussed in [i , by using the ray-shooting data structure [516111115] . 
with 0{n) time preprocessing, we can compute Vis(a) in 0(\Vis(a)\ logn) time, 
where \Vis(a)\ is the size of Vis(a). Note that it might be easier to compute 
|Vis(a)| by using both the ray-shooting data structure and our ray-rotating 
data structure in Lemma [1] 

The second issue is how to determine the next critical constraint of V that 
will be crossed by the point p as p moves from a to b. Suppose at the moment we 
know Vis{p) (initially, Vis(p) — Vis(a)). Let ft be the critical constraint that is 
crossed next by p. To determine ft, we first sketch an observation given in pQ. 

Denote by T(p) the shortest path tree rooted at p, which is the union of the 
shortest paths in V from p to all vertices of V. A vertex of V is in Vis(p) if and 
only if it is a child of p in T{p). For any child v of p in the tree T(p), define the 
principal child of v to be the vertex w among the children of v in T(p) such that 
the angle formed by the rays vu) and pv is the smallest among all such angles 
(see Fig. [5]). In other words, if we go from p to any child of v along the shortest 
path and we turn to the left (resp., right), then w is the first child of v that is 
hit by rotating counterclockwise (resp., clockwise) the line containing pv around 
v. 

To determine the next critical constraint ft, the following observation was 
shown in pQ. Two children of p in T(p) are consecutive if there is no other child 
of p between them in the cyclic order around p. 

Observation 1 pQ The next critical constraint ft is defined by two vertices of 
V that are either two consecutive children of p or one, say v, is a child of p and 
the other is the principal child of v. 

Based on Observation [I] Aronov et al. pQ maintained T(p) as p moves and 
used the balanced triangulation of V to determine the principal children. Their 
query algorithm takes <3(Hog 2 n) time, where k — \Vis(s)\, and the prepro- 
cessing time and space of their data structure are <3(n 2 log n) and 0(n 2 ), 
respectively. 

Here, we take a different approach, although we still use Observation [TJ 
Our data structure consists of the following: the ray-shooting data structure 
|5|6|11|15] . the ray- rotating data structure in Lemma [lj and a priority queue Q. 
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We assume Vis(a) has been computed. We use the ray-rotating data struc- 
ture to determine the principal children of all children of p, as follows. First of 
all, since we already know Vis(p) (initially p — a), we have all p's children in 
T(p), sorted cyclically around p. Note that we need not store the entire tree 
T(p). Consider any child v of p in T(p) (i.e., v is visible to p). Let w be the 
principal child of v that we are looking for. Consider the ray p(v) originating 
from v with the direction from p to v. By the definition of principle children, w 
is the first vertex of V visible to v that will be hit by the ray p{v) if we rotate 
p(v) around v along the direction that is consistent with the turning direction 
of the shortest paths from p to the children of v in T(p). It is easy to see that 
once we know the above rotation direction, we can obtain w in 0(log n) time by 
our ray-rotating data structure in Lemma [TJ 

To determine the above rotation direction, we only need to look at the re- 
lationship between the line containing p(v) and the two edges of V adjacent to 
v. Specifically, assume the line containing p(v) has the same direction as p{v). 
For example, if the two adjacent edges of v both lie to the left of this line (e.g., 
see Fig. [SJ, then we should rotate p(v) counterclockwise to determine w. The 
other cases can be determined in a similar manner. In summary, we can obtain 
the principle child of v in O(logn) time. Initially, p = a and we determine the 
principle children of all children of a in T(a) in 0(\Vis(a) \ logn) time since a 
has 0(\Vis(a)\) children in T(a). 

We use the priority queue Q to store the critical constraints specified in Ob- 
servation [T] that intersect the line segment s, where the key of each such critical 
constraint used in the priority queue Q is the position of its intersection with s. 
Initially when p — a, we compute the critical constraints defined by all pairs of 
consecutive children of p in T(p). Similarly, for each child v of T(p), we compute 
the critical constraint defined by v and its principal child. Note that the total 
number of these critical constraints is 0(\Vis(a)\). For each such critical con- 
straint, we check whether it intersects s, which can be done in 0(log n) time with 
the help of a ray-shooting query (we omit the details) . If the critical constraint 
intersects s, we insert it into Q; otherwise, we do nothing. Then, the first critical 
constraint in Q is the next critical constraint that p will cross as it moves. In 
general, after p crosses a critical constraint, p either sees one more vertex or sees 
one less vertex of V. In either case, there are only a constant number of insertion 
or deletion operations on Q. Specifically, consider the case when p sees one more 
vertex u (and an adjacent edge of u). By the implementation given in [2J, we can 
update the combinatorial representation of Vis{p) in constant time (i.e., insert 
u and the adjacent edge to the appropriate positions of the cyclic list of Vis(p)). 
After this, u becomes a child of p in the new tree T(p), and we can determine p's 
two other children, say, ui and 112, which are cyclic neighbors of u, in constant 
time. Then, for u\, we check whether the critical constraint defined by u and u\ 
intersects s, and if so, we insert it into Q. For U2, we do the same thing. Further, 
we compute the principal child of u in T(p), in O(logn) time, by the approach 
discussed above. For the other case where p sees one less vertex after it crosses 
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the critical constraint, we perform similar processing. After p arrives at the other 
endpoint b of s, we obtain the combinatorial representation of Vis(s). 

We claim that the above algorithm takes O(fclogn) time (with k — \Vis(s)\). 
Indeed, the initialization takes 0(\Vis(a)\logn) time. Clearly, |Vis(a)| = 0(k) 
since each vertex of V that is in Vis(a) also appears in Vis(s). If we consider 
every time when p crosses a critical constraint as an event, then each event takes 
O(logn) time. It has been shown in [TJ that the total number of events as p 
moves from a to & is O(k). Hence, the overall running time for computing Vis(s) 
is O(klogn). 

For the preprocessing, the ray-shooting data structure and the ray-rotating 
data structure both need 0(n) time and space to build. Further, in our query 
algorithm, the space used in the priority queue Q is always bounded by 0(k). 
We conclude this section with the following result. 

Theorem 1. For any simple polygon V , a data structure can be built in 0(n) 
time and 0(n) space, such that the visibility polygon Vis(s) can be computed in 
0(\Vis(s)\ log ro) time for any query line segment s in V. 

4 The Second Data Structure 

In general, the preprocessing of our second data structure is very similar to that 
in [2], and we make it faster by using better tools. Our improvement on the 
query algorithm is based on a number of new observations, e.g., a combinatorial 
bound of the "zone" of the line segment arrangements in simple polygons. For 
completeness, we first briefly discuss the approach in [2J. 

The preprocessing in [2J has several steps, whose running time is 0(n 3 logn) 
and is dominated by the first two steps. The other steps together take 0(n 3 ) 
time. We show below that the first two steps can be implemented in 0(n 3 ) time. 

The preprocessing in [2J first computes the visibility decomposition VD(V) 
of V. Although there may be fi{n 2 ) critical constraints in V, it has been shown 
[2] that any line segment in V can intersect only O(n) critical constraints, which 
implies that the size of VD(V) is 0(n 3 ) instead of 0(n A ). All critical constraints 
of V can be computed in 0(n 2 ) time, e.g., by the algorithm in [TO]. After 
that, to compute VD(V), we can use Chazelle and Edelsbrunner's algorithm 
[3] , which computes the planar subdivision induced by a set of m line segments 
in 0(m log m + /) time, where / is the number of intersections of the line seg- 
ments. In our problem, we have 0(n 2 ) critical constraints each of which is a line 
segment and the boundary of V has n edges. Therefore, by using the algorithm 
in [4], we can compute VD(V) in 0(n 3 ) time. Alternatively, an approach men- 
tioned in [15] can also be used to compute VD(V) in 0(n 3 ) time, and we omit 
the details. 

The second step of the preprocessing in [2] is to build a planar point location 
data structure on VD(V) in 0(n 3 logn) time. By the approaches in [OJ or [TO] , 
we can build such a point location data structure in 0(n 3 ) time. 

The remaining steps of our preprocessing algorithm are the same as those 
in [2], which together take 0(n 3 ) time. Hence, the total preprocessing time is 
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0(n 3 ). With the preprocessing, for each query point q in V , we can compute the 
visibility polygon Vis(q) of q in 0(\Vis(q) \ + logn) time. 

For a query segment s = ab, the query algorithm in [5] first computes Vis(a). 
Then, again, let a point p move on s from a to b. The algorithm maintains Vis(p) 
as p moves on s, initially with V^s(p) = Vis(a). Again, whenever p crosses a 
critical constraint, the combinatorial representation of Vis{j>) changes. Unlike 
our hrst data structure in Section [31 here we have VD(V) explicitly. Therefore, 
we can determine the next critical constraint in a much easier way. Specifically, 
the algorithm in [2] uses the following approach. Suppose p is currently in a cell 
of VD(V); then the next critical constraint crossed by p must be on the boundary 
of that cell. Since each cell is convex, we can determine this critical constraint in 
O(logn) time. The algorithm stops when p arrives at b. The total running time 
of the query algorithm is 0(/clogn), where k = \Vis(s)\. 

We propose a new and simpler query algorithm. We follow the previous query 
algorithmic scheme. The only difference is when we determine the next critical 
constraint that will be crossed by p, we simply check each edge on the boundary 
of the current cell that contains p, and the running time is linear in terms of 
the number of edges of the cell. Therefore, the total running time of finding 
all critical constraints crossed by p as it moves on s is proportional to the total 
number of edges on all faces of VD(V) that intersect s, and we denote by F(s) the 
set of such faces of VD(V). Let E(s) denote the set of edges of the faces in F(s). 
Then the total time of finding all critical constraints crossed by p is 0(\E(s)\). 
Note that the time of the overall query algorithm is the sum of the time for 
computing Vis(a) and the time for finding all critical constraints crossed by p. 
Since Vis(a) can be found in 0(|yis(a)| + logn) time, the running time of the 
query algorithm is 0(logn + |Uis(a)| + |i?(s)|). Recall that |Uis(a)| = 0(k). In 
Lemma [2] below, we will prove |-E(s)| = 0(k). Consequently, the query algorithm 
takes (3(logn + k) time and Theorem [3] below thus follows. 

Lemma 2. The size of the set E(s) is 0(k). 

Theorem 2. For any simple polygon V , we can build a data structure of size 
0(n 3 ) in 0(n 3 ) time that can compute Vis(s) in 0(\Vis(s)\ + logn) time for 
each query segment s in V . 

4.1 Proving Lemma [2] 

It remains to prove Lemma [21 Note that each edge of E{s) lies either on dV or 
on a critical constraint. We partition the set E(s) into two subsets Ei(s) and 
E 2 (s). For each edge of E(s), if it lies on dV, then it is in £a(s); otherwise, it is 
in E 2 (s). We will show that both = O(k) and \E 2 (s)\ = O(k) hold. 

Denote by C(s) the set of all critical constraints each of which contains at 
least one edge of E(s). 

Lemma 3. The size of the set C(s) is O(k). 
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Fig. 6. Illustrating the case where c inter- Fig. 7. Illustrating the case where c does 
sects s. not intersect s. 



Proof: Denote by V(s) the set of vertices of V visible to s. Clearly, |V(s)| < k. 
Consider an arbitrary critical constraint c £ C(s). To prove the lemma, we will 
charge c to a vertex of V(s). We will show that each vertex of V(s) will be 
charged at most a constant number of times, which will lead to the lemma. 

Assume the defining vertex of c is u and the anchor vertex of c is v. By the 
definition of C(s), c contains at least one edge of E(s). Depending on whether c 
intersects s, there are two cases. 

1. If c intersects s, then the defining vertex u of c must be visible to s (see 
Fig. [5]). To see this, let q be the intersection of c and s. Hence, q is visible to 
v. One may consider that the visibility between q and u is blocked by v. Due 
to our assumption that each endpoint of s is not collinear with two vertices 
of V, q is not an endpoint of s. Hence, there is always a point on s infinitely 
close to q that is visible to u, and thus u is visible to s. We charge c to its 
defining vertex u. 

2. If c does not intersect s (see Fig. [7]), then we show below that the anchor 
vertex v of c must be visible to s, and further, there are at most two critical 
constraints in C(s) such that they do not intersect s and their anchor vertices 
are v. We will charge c to v. 

We first prove that v is visible to s. Since c contains at least one edge of 
E(s), there must be a face / of VD(P) intersecting s and the boundary 
of / has an edge e contained in c. Let x be an arbitrary interior point of 
e and let y be an arbitrary point on s that is contained in / (see Fig. [7]). 
Since / is convex, xy is contained in /, i.e., x is visible to y and xy does 
not intersect any other critical constraint of V than c (at x). To prove y is 
visible to the vertex v, consider a point q on xy moving from x to y. We 
claim that v is always visible to q as q moves. Indeed, initially q is at x, and 
v is visible to x because x is on the critical constraint c and v is the anchor 
vertex of c. Suppose to the contrary v is not visible to q at some moment as 
q moves. Then, at some moment, vq must encounter a vertex of V, say, w. 
In other words, w is on vq. Then, the two vertices v and w define a critical 
constraint with v as the defining vertex and w as the anchor vertex, and the 
critical constraint intersects xy at q. Note that this critical constraint is not c 
because v is the anchor vertex of c. Hence, we obtain a contradiction because 
c is the only critical constraint that intersects xy. Therefore, we conclude 
that v is visible to y. 
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Fig. 8. Illustrating the critical constraints (the dotted line segments) with v as their 
anchor vertex that do not intersect s. 



Next, we prove that there are at most two critical constraints in C(s) such 
that they do not intersect s and their anchor vertices are v. Let C v denote 
the set of all critical constraints each of which has v as its anchor vertex and 
does not intersect s. Our goal is to prove that C v has at most two critical 
constraints in C(s). Note that each critical constraint in C v has v as an 
endpoint and its other endpoint is on dP. Hence, the critical constraints of 
C v partition V into \C V \ + 1 interior-disjoint regions and one region contains 
s (see Fig. [5]). Let R(s) be the region containing s. Clearly, R(s) has at most 
two critical constraints of C v , say c\ and C2 , on its boundary. We claim that 
for any critical constraint d 6 C v \ {ci,C2}, d cannot contain an edge of 
E(s). Indeed, assume to the contrary d contains an edge of E(s). Then, as 
discussed before, we can always find such two points x and y as in Fig. [7j 
Recall that xy is in V and xy does not intersect any other critical constraint 
of V than d . Since d is outside R(s), x G d is outside R(s). However, due to 
y S s and s C R{s), xy must intersect either c\ or C2, which contradicts with 
that xy does not intersect any other critical constraint of V than d . Hence, 
d cannot contain an edge of E(s) and d C(s). Therefore, we obtain that 
C v has at most two critical constraints in C(s). 

According to our discussion above, in the first case (i.e., c intersects s), we 
charge c to its defining vertex it, which is in V(s). In the second case (i.e., c does 
not intersect s), we charge c to its anchor vertex v, which is also in V(s). An 
observation in [2] shows that for any line segment in V, for any vertex u of V, the 
line segment intersects at most two critical constraints with u as their defining 
vertex. Therefore, for any vertex u of T 3 , u can be charged at most twice as a 
defining vertex. On the other hand, we have shown that, as an anchor vertex, v 
has at most two critical constraints in C(s) that do not intersect c. Therefore, 
for any vertex v of V, v can be charged at most twice as an anchor vertex. Hence, 
any vertex in V(s) can be charged at most four times, twice as an anchor vertex 
and twice as a defining vertex. In other words, |C(s)| < 4 • |V(s)| < 4 • k. 

The lemma thus follows. □ 
In the next lemma, we bound the size of the subset E\ (s) . 

Lemma 4. The size of the set Ei(s) is 0(k). 

Proof: Denote by V(s) the set of vertices of V visible to s. Clearly, |V(s)| < k. 
Consider an edge e in Ei(s). To prove the lemma, we will charge e either to a 
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vertex of V(s) or to a critical constraint of C(s). We will also show that each 
vertex of V(s) will be charged at most twice and each critical constraint of 
C(s) will be charged at most four times. Consequently, due to |V(s)| < k and 
|C(s)| = 0(k) (by Lemma[3]), the lemma follows. 

By the definition of Ei(s), e is on an edge of P. If e has an endpoint that is 
a vertex of V, say u, then clearly u is visible to s. We charge e to u. Otherwise, 
both endpoints of e are endpoints of some critical constraints, and we charge e 
to an arbitrary one of the two such critical constraints. 

For each vertex of V, it has two adjacent edges in V, and therefore, it has at 
most two adjacent edges in E±(s). Hence, each vertex of V(s) can be charged at 
most twice. On the other hand, each critical constraint has two endpoints, and 
each endpoint is adjacent to at most two edges in Ei(s). Therefore, each critical 
constraint of C(s) can be charged at most four times. □ 

To prove LemmaHJ it remains to show |i?2(s)| = 0(k). To this end, we discuss 
a more general problem, in the following. 

Assume we have a set S of line segments in V such that the endpoints of 
each such segment are on dV . Let A be the arrangement formed by the line 
segments of S and the edges of dV. For any line segment s in V (the endpoints 
of s need not be on dV), the zone of s is defined to be the set of all faces of A 
that s intersects. Denote by Z(s) the zone of s. For each edge of a face in A, it 
either lies on a line segment of S or lies on dV] if it is the former case, we call 
the edge an S-edge. We define the complexity of Z(s) as the number of 5-edges 
of the faces in Z(s) (namely, the edges on dV are not considered), denoted by 
A. Our goal is to find a good upper bound for A. By using the zone theorem for 
the general line segment arrangement [5], we can easily obtain .A = 0(|S , |a(|S'|)), 
where a(-) is the functional inverse of Ackermann's function |13j . 

Denote by S s the set of line segments in S that intersect Z(s), i.e., each 
segment in S s contains at least one S'-edge of Z(s). Let m = \S S \ (note that 
m < 15*1). By using the property that each segment in S has both endpoints on 
dV, we show A = 0(m) in Theorem [3] below, which we call the zone theorem. 
The proof of Theorem [3] is given in Section 14.21 

Theorem 3. The complexity of Z(s) is 0(m). 

Now consider our original problem of proving |i?2(s)| = 0(k). By using the 
zone theorem, we have the following corollary. 

Corollary 1. The size of the set i?2(s) is 0(k). 

Proof: The set i?2(s) consists of all edges of E(s) that lie on the critical con- 
straints. Recall that each critical constraint is a line segment in V with both 
endpoints on dV. Consider the arrangement formed by all critical constraints 
of V and dV. The complexity of the zone Z(s) of the query segment s in this 
arrangement is exactly |i?2(s)|- Let C'(s) be the set of critical constraints of V 
each of which contains at least one edge in E%{s). Then, by the zone theorem 
(Theorem El), we have |^ 2 (s)| = 0(\C'(s)\). Note that C'{s) C C(s). Due to 
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\C(s)\ = 0{k) (Lemma[3]), we have \E 2 {s)\ = 0(k). The corollary thus follows. 
□ 

Lemma 2] and Corollary [T] together lead to Lemma U 
4.2 Proving the Zone Theorem (i.e., Theorem [3]) 

This subsection is devoted entirely to proving the zone theorem, i.e., Theorem 
[3] All notations here are the same as defined before. 

We partition the set S s into two subsets: Si and . For each segment in S s , 
if it does not intersect the interior of s, then it is in Si; otherwise, it is in Sg. 
Let mi = \Sl\ and rri2 = \S^\. Hence, m = m\ + 7712. Consider the arrangement 
formed by the line segments in Si and dV . Since no segment in Si intersects 
the interior of s, s must be contained in a single face of the above arrangement 
and we denote by F s that face. For each edge of F s , if it lies on a segment of S, 
we also call it an S-edge. Note that the edges of F s that are not S-edges are all 
on dV. 

Lemma 5. The number of S-edges of the face F s is 0{mi); the shortest path 
in V between any two points in F s is contained in F s . 

Proof: For each segment s' in Sj, since both endpoints of s' are on dV, s' 
partitions V into two simple polygons and one of them contains s, which we 
denote by V(s'). It is easy to see that the face F s is the common intersection 
of V(s')'s for all s' in Si- To prove the lemma, it is sufficient to show that each 
segment s' in Si has at most one (maximal) continuous portion on the boundary 
of F s , as follows. 

For any two points p and q'mV, denote by n(p, q) the shortest path between 
p and q in V . Note that since V is a simple polygon, 7r(p, q) is unique. We claim 
that for any two points p and q in the face F s , ir(p, q) is contained in F s . Indeed, 
suppose to the contrary ir(p,q) is not contained in F s . Then, 7r(p, q) must cross 
the boundary of F s . Since ir(p,q) cannot cross the boundary of V , ir{p,q) must 
cross an S'-edge of F s , and we assume s' is the segment in Si that contains such 
an S'-edge. This implies that n(p,q) is also not contained in the polygon V(s'). 
Recall that the line segment s' partitions V into two simple polygons and one of 
them is V(s'). It is easy to show that for any two points in V(s'), their shortest 
path in V must be contained in V(s'). Therefore, we obtain a contradiction. 
Hence, our above claim is true. 

Now assume to the contrary that a segment s' in Si has two disjoint maximal 
continuous portions on the boundary of F s . Let p and q be two points on these 
two portions of s', respectively. Thus, both p and q are in F s . Since these are 
two discontinuous portions of s' on the boundary of F s , the line segment pq 
is not contained in F s . Since pq is on s', the shortest path Tr(p,q) is pq. But 
this means 7r(p, q) is not contained in F s , which incurs a contradiction with our 
previous claim that n(p, q) must be contained in F s . Hence, we obtain that each 
segment s' in Si has at most one continuous portion on the boundary of F s , and 
consequently, the number of S-edges of the face F s is 0(mi). □ 

Lemma |6] below shows a property of the face F s . 
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Lemma 6. For any line segment s' in V with both endpoints on dV, s' has at 
most one (maximal) continuous portion intersecting F s ; consequently, s' inter- 
sects the interior of at most two edges of F s . 

Proof: Assume to the contrary that s' has two disjoint maximal continuous 
portions intersecting F s . Let p and q be two points on these two portions of s', 
respectively. Thus, both p and q wee in F s . Clearly, the line segment pq is not 
contained in F s . Since pq is on s' , pq is the shortest path ir(p, q) between p and q 
in V . But this means 7r(p, q) is not contained in F s , which incurs a contradiction 
with Lemma [51 Hence, the lemma holds. □ 
For each 5-edge of Z{s), it lies either on a segment in S). or on a segment in 
Sf; we call it an S^-edge if it lies on a segment in S\ and an S^-edge otherwise. 
Due to m = mi + m 2 , our zone theorem is an immediate consequence of Lemma 
[7] below. Note that we can obtain the zone Z(s) of s by adding the segments of 
Sg to F s . To prove Lemma [71 we use induction on wi2, i.e., \S^\. The approach 
is very similar to that in [7] used for line arrangements. Here, although we have 
line segments instead of lines, the property that each line segment has both 
endpoints on dV makes the approach in [7J applicable with some modifications. 

Lemma 7. There are 0(7712) S^-edges and 0(m\ + 777,2) Sl-edges in the zone 
Z(s). 

Proof: Without loss of generality, assume the segment s is horizontal. It is easy 
to see that each S^-edge bounds one face of Z(s) and each S^-edge bounds two 
faces of Z(s) (one lies on its right and the other lies on its left). For each Sl- 
edge, we say it is a left bounding S 2 s -edge for the face lying on the right of it 
and a right bounding S 2 -edge for the face lying on the left of it. Below we will 
prove that the number of left bounding S^-edges of the faces in Z(s) is 0(7772). 
Analogously, the number of right bounding S^-edges of the faces in Z{s) is also 
0(7772). In addition, we will also show that the number of S^-edges of Z{s) is 
0(?77i + 7772) = 0(777). 

Our proof is by induction on 7772. Consider the base case with 7772 = 1. 
Denote by s' the only segment in S 2 . Note that the face F s has no S 2 -edges 
on its boundary and has 0(tt7i) S^-edges by Lemma [5] In light of Lemma [51 
s' has at most one maximal continuous portion intersecting the face F s and s' 
intersects the interior of at most two 5* -edges of F s . Therefore, after we add s' 
to F s , the number of S^-edges of Z(s) increases by at most two and the number 
of left bounding S^-edges increases by at most one. 

Consider the general case of 7772 > 1. Let s' be the segment in S 2 that 
intersects s at the rightmost position among all segments in S 2 . We first consider 
the case when this segment s' is unique. By induction, the zone of s has c- (7772 — 1) 
left bounding S^-edges and 0(7771+7772 — 1) S^-edges, for some constant c, without 
considering the segment s' . Now consider adding s' . First, by Lemma [51 s' has at 
most one maximal continuous portion intersecting the face F s and s' intersects 
the interior of at most two S^-edges of the zone Z(s); therefore, the number of 
Sl-edges increases by at most 2. Second, the number of left bounding S^-edges 
increases in two ways: there are new left bounding S^-edges on s' and there are 
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Fig. 9. The shaded region is R(v), which is not in the zone of s. 



old left bounding £3 -edges that are split by s'. Let v be the first intersection 
point of s' with another segment in S% above s, and let w be the first intersection 
point of s' with another segment in Sg below s (e.g., see Fig.[9]). We assume both 
v and w exist since otherwise the analysis is even simpler. The segment vw on 
s' becomes a new left bounding S^-edge. In addition, s' splits a left bounding 
S^-edge at v and at w, respectively. Hence, the number of the left bounding 
S^-edges increases by three. We claim that there is no other increase for the 
number of left bounding S^-edges. 

Indeed, consider the part of s' above v. Let s" be the segment in Sg that 
intersects s' at v. Let p' be the endpoint of s' above v and p" be the endpoint 
of a" above v. Note that both p' and p" are on dV. Consider the region R(v) 
above v enclosed by vp' , vp" , and the portion of dV between p' and p" (e.g., see 
Fig. in]). Clearly, the region R(v) is not in the zone of s. Further, R(v) lies on the 
right of vp', and thus vp' cannot contribute any left bounding 5^-edges to Z(s). 
In addition, if a left bounding S^-edge e that was in the zone Z(s) (before s' is 
added) is intersected by s' somewhere above v, then the part of e to the right of 
s' (i.e., the part of e in the region R(v)) is not in the zone Z(s) any more after 
s' is added. Hence, there is no increase in the number of left bounding -edges 
due to such an intersection. 

In a similar way, we can show that the part of s below w does not increase the 
number of left bounding 5;:-cdges in the zone Z(s). Therefore, after s' is added, 
the total increase of the number of left bounding ^-edges is at most three. 

We discuss above the case when s' is the only segment in through the 
rightmost intersection on s. If there is more than one such segment, then we 
take an arbitrary such segment as s'. By a similar analysis as that above and 
that in [7] , we can show that the total increase in the number of left bounding 
S^-edges is at most five. We omit the details. 

We conclude that there are 0(m2) .S^-edges and 0(m\ +1112) Sg-edges in the 
zone Z(s). The lemma thus follows. □ 

5 Conclusions 

In this paper, we propose two new data structures for the weak visibility query 
problem on a simple polygon, which improve upon the previous work |l|2|3j . 



15 



Some results (e.g., the ray- rotating data structure and the zone theorem) may 
be of independent interest. 

For the O(fclogn) time queries, our first data structure is clearly optimal. 
For the 0(k + logn) time query, however, an open question is whether a data 
structure of sub-cubic preprocessing time and space is possible. 
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